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>• ' Abstract 

' The experimental reconstruction of the Higgs self-energy potential is essential to a 

. verification of the Higgs boson's role in spontaneous electroweak symmetry breaking. 

^) • The first step towards this goal, the measurement of the triple Higgs self-coupling, 

o ■ 

' can possibly be accomplished at the next generation of linear colliders. Here we 

Q^l present a background study of the most promising channel, double Higgs-strahlung 

Qh! off a Z boson, e+e" — > HHZ, with the subsequent decay H bb, and evaluate the 

(D . 

, feasibility of its measurement. 

> 

H : 1 Introduction 



To verify whether or not the Higgs mechanism is responsible for spontaneous electroweak 
symmetry breaking as expected by the standard model (Sm), or indeed, the minimal su- 
persymmetric standard model (Mssm), one must perform an experimental reconstruction 
of the Higgs self-energy potential. This reconstruction requires the measurement of the 
Higgs boson mass(es) and the triple and quartic Higgs self-couplings. The triple Higgs 
self-coupling may become accessible to direct measurement at an electron-positron linear 
collider (Lc) with centre-of-mass energy at the TeV scale. In this study we compare 
the signal of the most promising channel with its dominant electroweak (Ew) and QCD 
backgrounds to assess the feasibility of its measurement. 



^Talk given at the Fourth Workshop of the "2nd Joint ECFA/DESY Study on Physics and Detectors 
for a Linear Electron- Positron Collider," Oxford, UK, 20-23 March 1999. 
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We restrict ourselves to a discussion of the triple Higgs self-coupling in the Sm0 and in 
particular double Higgs-strahlung off Z bosons, in the process e^e^ — > HHZ . We adopt 
the H ^ bb decay channel over the Higgs mass range Mh ^ 140 GeV and assume very 
efficient tagging and high-purity sampling of b quarks. Then the backgrounds to the Xrhh 
measurement are primarily the 'irreducible' Ew and QCD backgrounds e~^e~ bbbbZ. 

The Ew background is of 0{al^) away from resonances, but can, in principle, be 
problematic due to the presence of both Z vectors and Higgs scalars yielding bb pairs. In 
contrast, the QcD background is of C(ag^a^). Here, although there are no heavy objects 
decaying to bb pairs, the production rate itself could give difficulties due to the presence of 
the strong coupling. In addition, the double Higgs-strahlung process (see Fig. ^ contains 
diagrams proceeding via an HHZ intermediate state but not dependent on Xhhh, as well 
as a diagram sensitive to the triple Higgs self-coupling (the right-hand graph of Fig. §). 




Figure 1: Diagrams contributing at lowest-order to e+e^ — > &(3)fe(4)&(5)&(6)^(7) via purely 
Ew interactions containing two Higgs bosons in the intermediate state. Note that only 
the far right diagram contains a dependence on Xhhh- 

The plan of the paper is as follows. The next section details the procedure adopted 
in computing the relevant scattering amplitudes. Section |^ displays our numerical results 
and contains our discussion. Finally, in the last section, we summarize and conclude. 

2 The calculation 

The signal process is rather straightforward to calculate in the case of on-shell HHZ pro- 
duction (see Refs. 0, ^ for analytic expressions of the matrix elements (Mes) and further 
discussions of related signal processes). The Ew background is more complex, deriving 

^Some of the trilinear couplings of the MSSM may prove to be accessible even at the Lhc via resonant 
decay of a heavy Higgs, eg. H — > hh, and are under investigation elsewhere [Q. 
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from many graphs with different resonant structures. In order to perform an efficient in- 
tegration we have grouped the Feynman diagrams into different collections of diagrams 
with identical (non-)resonant structure. This categorization allows one to compute each 
of the topologies separately, with the appropriate mapping of variables, thus optimizing 
the accuracy of the numerical integration. Furthermore, one is able to assess the relative 
weight of the various subprocesses in the full scattering amplitude, giving added insight 
into the fundamental dynamics. In contrast the QcD background contains fewer diagrams, 
with only five different (non-) resonant topologies. This makes integration much simpler 
than in the Ew case and it can be carried out with percent accuracy directly over the full 
Me using standard multichannel Monte Carlo methods. For further details and numerical 
inputs see Ref. 

We assume that total and differential rates are those at the partonic level, as we identify 
jets with the partons from which they originate. To resolve the final state b (anti)quarks 
as separate systems, we impose the following acceptance cuts: E{b) > 10 GeV on the 
energy of each b (anti)quark and cos(6, b) < 0.95 on the relative separation of all possible 
26 combinations. We further assume that b jets are distinguishable from light-quark and 
gluon jets (by using, for example, /i- vertex tagging techniques). However, no efficiency to 
tag the four b quarks is included in our results. The Z boson is treated as on-shell and no 
branching ratio is applied to quantify its possible decays. In practice, in order to simplify 
the treatment of the final state, one may assume that the Z boson decays leptonically (that 
is, Z —>■ , with i = e, fi, r) or hadronically into light quark jets (that is, Z qq, with 
q^ b). Also, we have not included Initial State Radiation (Isr) in our calculations. 

3 Results 

The total signal cross section is plotted as a function of in the top-left frame of 
Fig. ^, for three centre-of-mass (Cm) energies. Even at low Higgs masses where both the 
production and decay rates are largest, the signal is rather small. In fact, the signal is below 
0.2 femtobarns for all energies from 500 to 1500 GeV, although this can be doubled simply 
by polarizing the incoming electron and positron beams. Thus, as already recognized 
in Ref. 0, where on-shell production studies of the signal were performed, luminosities 
of the order of one inverse attobarn need to be collected before statistically significant 
measurements of Xhhh can be performed. 
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The decrease of the signal with increasing Higgs mass is due to suppression of the 
H ^ bb decay channel. Above Mh ~ 140 GeV it becomes overwhelmed by the opening 
of the off-shell H decay (see, for example, Fig. 1 of Ref. 0). In contrast, 

the production cross section for e+e~ HHZ without specifying the subsequent decay 
is much less sensitive to Mfj 0. In addition, because the signal is an annihilation process 
proportional to a larger Cm energy (-Ecm) tends to deplete the production rates, as 
long as Ectn > 2Mh + Mz- When this is no longer true, e.g., at 500 GeV and Mh ^ 140 
GeV, phase space suppression can overturn the 1/s propagator effects. This is evidenced 
by the crossing of the curves for 500 and 1000 GeV in Fig. ^. 





Figure 2: Top-left: cross sections in femtobarns for the signal at three different collider 
energies: 500, 1000 and 1500 GeV. Top-right (Bottom-left) [Bottom-right]: cross sections in 
femtobarns for the signal versus the Ew and QCD backgrounds at 500(1000) [1500] GeV. 
Our acceptance cuts in energy and separation of the four b quarks have been implemented. 

Fig. also shows the background processes plotted with respect to Mh for the three 
Cm energies. As anticipated, the Ew background is problematic due to its resonance 
structures, whereas the QcD background presents no such difficulty and does not signif- 
icantly obscure the signal, despite the strong coupling constant. Our categorization of 
the Ew background into different resonant topologies now facilitates a closer examination. 
In particular we observe that only four sub-processes dominate the background. Generic 



4 



Feynman diagrams of these sub-processes can be seen in Fig. together with their con- 
tribution to the total Ew background rate, as a function of the Higgs boson mass. All 
other Ew sub-processes are much smaller, rarely exceeding 10~^ femtobarns, and have 
little effect on the kinematics or magnitude of the background. 




sub-process 3 sub-process 4 i,^ (c^v) 



Figure 3: On the left, generic Feynman diagrams representing the four dominant Ew 
topologies. Note that only sub-process 1 does not contain a Higgs boson in the intermediate 
state. The different contributions of these sub-processes to the total Ew background are 
displayed in the plot to the right, as a function of the Higgs boson mass. 



The QcD background is dominated by e"'"e^ —>■ ZZ production with one of the two 
Z bosons decaying hadronically into four b jets. Also significant is single Higgs-strahlung 
production (off a Z) with the Higgs scalar subsequently decaying into bbbb via an off-shell 
gluon. The contributions of the other diagrams, which do not resonate, are typically one 
order of magnitude smaller than the ZZ and ZH mediated graphs, with the interferences 
smaller still (and generally negative). 

We now investigate several differential spectra, to find kinematic cuts which will sup- 
press the backgrounds. The distributions in E{b) and cos(6, b) cannot be further ex- 
ploited after the acceptance cuts are made. Instead we consider the invariant masses of b 
(anti)quark systems: for 2b systems where the 6 jets come from the same production vertex 
('right' pairing) or otherwise ('wrong' pairing); for 36 systems, and for the 46 system. The 
spectra for the 2b and 46 systems are shown in Fig. ^ (the 36 spectra are less useful and are 
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not shown here). It is clear that the narrow Higgs peak[| in the M^lbb) distribution, can 
be exploited to reduce the QcD background. Although the Ew process also displays a res- 
onance at Mh, only one 2b invariant mass will peak here compared to two combinations in 
the signal. The Ew background can therefore also be cut away. Finally, one may require 
that none of the 26 invariant masses reproduce a Z boson, provided that the invariant 
mass resolution of di-jet systems is at least as good as the difference {Mh — Mz)/2, in 
order to resolve the Z and H peaks. The M{bbbb) spectra is also useful in distinguishing 
signal from background since clearly M{bbbb) must always be greater than 2Mh for the 
signal process, while it can be lower for both background processes (especially for the QcD 
background) . 
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Figure 4: Differential distributions in invariant mass of multi-jet systems containing two 
or four h (anti) quarks. Also shown in the bottom right plot is the (cosine of) the angle of 
the 46 system with respect to the beam axis. The Cm energy is 500 GeV and the Higgs 
mass 110 GeV. Our acceptance cuts in energy and separation of the four h quarks have 
been implemented. 



In addition to the different resonance structures of signal and background one can 
exploit the dominantly t-channel nature of the backgrounds as compared to the 5-channel 



■^Recall that for Mh = 110 GeV one has Th ~ 3 MeV. The Higgs resonances in Fig. ^ are smeared by 
incorporating a 5 GeV bin width. 
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signal. In Fig. ^ we also show the cosine of the polar angle (i.e. with respect to the beam 
axes) of the four b quark system (or, indeed, the real Z). Notice that the backgrounds are 
much more forward peaked than the signal. The QCD events are predominantly e~^e~ —>■ 
ZZ production followed by the decay of one of the gauge bosons into four b quarks. The 
ZZ pair is produced via t, u-channel graphs so the four b quarks are preferentially directed 
forwards and backwards into the detector. In contrast, the signal is entirely s-channel 
resulting in more centrally produced b jets. The EW background has a more complicated 
structure but is still sizably dominated by forward production. A similar effect is seen for 
cos{bbb) and cos{bb), allowing one to separate signal and backgrounds events efficiently. 

The transverse momentum distributions of the 46, 36 and 2b systems were also exam- 
ined. However the distributions for the signal and the Ew background proved too similar 
to allow their use as efficient kinematic variables for cuts. 

The main features of distributions studied above are rather stable to variation of the 
Cm energy or Higgs mass (within the ranges under discussion). We can therefore optimize 
the S/ B ratio by imposing the cuts: 

\M{bb) — Mh\ < 5 GeV (on exactly two combinations of 26 systems), 

|M(66) — Mz\ > 5 GeV (for all combinations of 26 systems) 

M(6666) > 2Mh, | cos(26, 36, 46)| < 0.75. (1) 

In enforcing these constraints, we assume no 6 jet charge determination. 

The signal and background cross-sections after the implementation of the selection cuts 
can be seen in Fig. ^. Both background rates are greatly reduced while a large portion of the 
original signal is maintained. This results in a S/B ratio which is enormously large for not 
too heavy Higgs masses. For example, at Centre-of-mass energies Ecm = 500(1000) [1500] 
GeV and for Mh = 110 GeV, one finds S/B = 25(60) [104]. This remarkable suppression 
of the backgrounds comes largely from the invariant mass cuts on Mhh. In fact, they are 
crucial not only in selecting the Mh resonance of the signal, but also in minimizing the 
signal rejection around Mz when mispairings occur: notice the shoulder at 90 GeV of the 
Mw{bb) signal spectrum. 

A number of caveats to our analysis apply. Firstly, the value we have adopted for 
the resolution is rather high, considering the large uncertainties normally associated with 
the experimental determination of jet angles and energies, though not unrealistic in view 
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Figure 5: Top-left: cross sections in femtobarns for the signal at three different collider 
energies: 500, 1000 and 1500 GeV. Top-right (Bottom-left) [Bottom-right]: cross sections in 
femtobarns for the signal versus the Ew and QcD backgrounds at 500(1000) [1500] GeV. 
Acceptance and selection cuts have been implemented. 



of the most recent studies [0|. The ability of the actual detectors in guaranteeing the 
performances foreseen at present is thus crucial for the feasibility of dedicated studies of 
double Higgs-strahlung events at the LC. Furthermore, one must consider the efficiency of 
tagging the b quarks necessarily present in the final state, particularly in the case in which 
the Z boson decays hadronically. Given the high production rate of six jet events from QcD 



IP and multiple gauge boson resonances |]T0[ in light quark and gluon jets, it is desirable 
to resort to heavy fiavour identification in hadronic samples. However, the poor statistics 
of the HHZ signal requires a judicious approach in order not to deplete it below detection 



level. According to recent studies [|TT|, the two instances can be combined successfully, as 
efficiencies for tagging bb pairs produced in Higgs decays were computed to be as large as 
e^fe ~ 90%, with mis-identification probabihties of light (charmed) quarks as low as egg(cc) ~ 
0.3(4)% (and negligible for gluons). If such a projection for the Lc detectors proves to be 
true, then even the requirement of tagging exactly four b quarks in double Higgs-strahlung 
events might be statistically feasible, thus suppressing the reducible backgrounds to really 
marginal levels [jl2| . One should also bear in mind that experimental considerations, such 
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as the performances of detectors, the fragmentation/hadronization dynamics and a reahstic 
treatment of the Z boson decays, are also important when determining what cuts should 
be made. Such considerations are beyond the scope of this paper, and are under study 
elsewhere |12[ . 



Finally, the number of signal and backgrounds events seen per inverse attobarn of 
luminosity at i^cm = 500, 1000, and 1500 GeV, with Mn = 110 GeV, can be seen in 
Tab. ^ One could relax one or more of the constraints we have adopted to try to improve 
the signal rates without letting the backgrounds become unmanageably large. However, 
such a relaxation only marginally effects the signal while significantly reducing the S/B 
ratio, and should only be done if high luminosities cannot be obtained. Kinematic fits can 
also help in improving the S/ B ratio [O . 





Number of Events per ah ^ after selection cuts 


Ecm = 500 GeV 


= 1000 GeV 


= 1500 GeV 


signal 


26 


40 


34 


Electroweak 


1.0 


0.6 


0.3 


QCD 


0.032 


0.026 


0.016 



Table 1: The number of events for signal and backgrounds per inverse attobarn of lumi- 
nosity after selection cuts for centre-of-mass energies of 500, 1000 and 1500 GeV, a Higgs 
mass of 110 GeV, and with polarized electron and positron beams. 



4 Summary 

In conclusion, the overwhelming irreducible background from Ew and QcD processes of 
the type e"'"e~ bbbbZ to double Higgs production in association with Z bosons and decay 
in the channel H — ^ bb, i.e., e'^e^ HHZ bbbbZ, should easily be suppressed down to 
manageable levels by simple kinematics cuts: e.g. in invariant masses and polar angles. 

The number of signal events is generally rather low, but will be observable at the Lc 
provided that it provides very high luminosity, excellent b tagging performances, and high 
di-jet resolution. As advocated in Ref. [|12|, one also requires a good forward acceptance 



for jets since single jet directions in the double Higgs-strahlung process can stretch up to 
about 20° in polar angle. 
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